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Based on multiyear INTEGRAL observations of SS433, a composite IBIS/ISGRI 18-60 keV 
orbital light curve is constructed around zero precessional phase = 0, which corresponds to 
a maximum separation of the moving emission lines originated in sub-relativistic jets from the 
source. It shows a peculiar shape characterized by a significant excess near the orbital phase 
0«ri = 0.25, which is not seen in the softer 2-10 keV energy band. Such a shape is likely to be 
due to a complex asymmetric structure of the funnel in a supercritical accretion disk in SS433. 
The orbital light curve at 40-60 keV demonstrates two almost equal bumps at phases ^ 0.25 
and ~ 0.75, most likely due to nutation effects of the accretion disk. The change of the off- 
eclipse 18-60 keV X-ray flux with the precessional phase shows a double-wave form with strong 
primary maximum at i/^,- = and weak but significant secondary maximum at = 0.6. A 
weak variability of the 18-60 keV flux in the middle of the orbital eclipse correlated with the 
disk precessional phase is also observed. The joint analysis of the broadband (18-60 keV) orbital 
and precessional light curves obtained by INTEGRAL confirms the presence of a hot extended 
corona in the central parts of the supercritical accretion disk and constrain the binary mass ratio in 
SS433 in the range 0.5 > ^ > 0.3, confirming the black hole nature of the compact object. Orbital 
and precessional light curves in the hardest X-ray band 40-60 keV, which is free from emission 
from thermal X-ray jets, are also best fitted by the same geometrical model with hot extended 
corona at q ^ 0.3, stressing the conclusions of the modeling of the broad-band X-ray orbital and 
precessional light curves. 
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1. Introduction 



SS433 is a unique galactic steadily superaccreting microquasar with mildly relativistic (v = 
0.26c), precessing jets located at a distance of 5.5 kpc 0, ^ ||, 0, ||]. The system exhibits three pho- 
tometric and spectral periodicities related to precession (Ppiec = 162'^'. 5), orbital (/"orb = 13'^.082) 
and nutation (P^ut = 6^.28) periods [pi]]. Despite of wealth of observations, the nature of the com- 
pact star in SS433 remains inconclusive. The presence of absorption lines in the optical spectrum of 
the companion [10, |l2|] suggests its spectral classification as A7Ib supergiant. Assuming these 
lines to be produced in the optical star photosphere, their observed orbital Doppler shifts would 
correspond to the mass ratio of compact (Mv) and optical (My) star q = M^/My ~ 0.3 ±0.11 and 
masses A/^ = (4.3 ± 0.%)Mq, My = (12.3 ± 3.3)M0, respectively, pointing to the black hole nature 
of the compact star. 

Modeling of all INTEGRAL eclipses of the source available before 2010 [Q] using a purely 
geometrical model yeilded independent constraints on the binary mass ratio q = 0.25 — 0.5 with the 
most probable value q = 0.3, suggesting the mass of the compact companion ^ 5.3Mq and the 
optical star My ~ 17. VM© for the observed optical star mass function fy = O.268M0. This places 
SS433 among black-hole high-mass X-ray binaries. Thus SS433 can be the only known example of 
galactic massive X-ray binary at an advanced evolutionary stage [^ ^ with supercritical accretion 
[ jl6| ] onto a black hole [^. Its study in different spectral bands provides invaluable information for 
theory of evolution of binary star and the formation of relativistic jets. 

The basic picture of hard X-ray emitting regions, as emerged from analysis of X-ray data 
[^, ^, ^, ^, 13], includes hot X-ray jet propagating through a funnel in the supercritical accretion 
disk, filled with hot scattering medium (a corona). The X-ray spectrum of SS433 in the 3-100 keV 
range can be fitted by two-component model (thermal X-ray emission from the jet and thermal 
comptonization spectrum from corona) elaborated in [13|. The scattering corona parameters are: 
Tcor — 20 keV, Thomson optical depth Tj- ~ 0.2 and mass outflow rate in the jet Mj = 3 x 10^^ g/s. 
This parameters suggest the coronal electron number density around 5 x lO'^ cm^^, which is typ- 
ical in the wind outflowing with a velocity of v ~ 3000 km/s from a supercritical accretion disk 
with mass accretion rate onto the compact star M 10^^ M0/yr at distances 10^^ cm from the 
center, where a Compton-thick photosphere is formed [|]|. The size of the disk photosphere was 
independently estimated from measurements of fast optical aperiodic variability [jl]]. 

Here we analyze hard X-ray eclipses of SS433 near the T3 moment in combination with the 
precessional variability as observed by INTEGRAL, and interpret them in terms of our multicom- 
ponent geometrical model (see [^ for more detail). 



2. Composite X-ray light curve and jet nutation effect 

Previous dedicated INTEGRAL observations of hard X-ray eclipse are summarized in Table 
|TJ They were all concentrated around precessional phase zero ('the T3 moment' in terms of the 
kinematic model of SS433 [p^), where the accretion disk is maximum opened to the observer and 
the X-ray flux from the source is the highest. 

The composite IBIS/ISGRI 1 8-60 keV light curve of the primary X-ray eclipse at precessional 
phase zero is shown in Fig. |l] (panel 1). Data were analyzed using the IKI INTEGRAL data 
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Table 1: Dedicated observations of SS433 by INTEGRAL 



Set INTEGRAL orbits Dates Precessional phase xj/pr 



I 


67-70 


May 2003 


0.001-0.06 


II 


555-556 


May 2007 


0.98-0.014 


III 


608-609 


October 2007 


0.956-0.99 


IV 


612-613 


October 2007 


0.030-0.064 


V 


722-723 


September 2008 


0.057-0.091 


VI 


984 


November 2010 


0.87-0.89 




987 


November 2010 


0.93-0.94 


VII 


1040-1041 


April 2011 


0.91-0.95 




Figure 1: From left to right: 1) Composite IBIS/ISGRI 18-60 keV X-ray eclipse light curve around zero 
precessional phase Xj/^O. 2) Binned X-ray light curve (A0 = 0.02) with phase intervals for spectral analysis. 
3) Phase-resolved IBIS/ISGRI spectra of SS433 within chosen orbital phase intervals, l-a errors are shown. 



processing code described in [[Isll. The eclipse light curve averaged within orbital phase bins 
= 0.02 is shown in Fig. [T] (panel 2). 3a-flux errors are indicated. The 18-40 and 40-60 keV 
orbital and precessional light curves with the corresponding hardness ratios are shown in Fig. || 
and Fig. ||, respectively. 

A very significant flux excess at the orbital phase (p 0.25 is observed in the composite 18-60 
light curve (Fig. |l|, panel 2; Fig. ^, upper light curve) after the eclipse relative to the phase 0.75 
before the eclipse. However, on the 40-60 keV light curve two maxima with similar amplitude 
are clearly seen at both orbital phases 0.25 and 0.75. These ~ 10% sine-like variability at twice 
the orbital period superimposed on the orbital light curve is mostly likely to be due to the jet 
nutation. The tidal nutation occurs twice the synodic period (6.28 d in the case of SS433). The 
effect must be maximal at the zero precessional phase when the binary system is observed in 
quadratures (i.e. at the binary phases 0.25 and 0.75). Due to nutation the jet changes its inclination 
to the line of sight by 6 degrees, suggesting the change in the projection area of the jet funnel 
AS'/S' ~ Aij/ij = 0.1 (here ij 59" is the jet inclinatin angle at the T3 precessional phase). Thus, 
the composite INTEGRAL light curve allowed us to see for the first time the jet nutation in the 
40-60 keV band. 

2. Orbital phase-resolved spectroscopy. The accumulated data allowed us to make for the first 
time the orbital phase-resolved X-ray spectroscopy. Five orbital phase intervals chosen for spectral 
analysis are shown by the vertical dashed lines in Fig. |T] (middle panel). The obtained X-ray 
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Figure 2: IBIS/ISGRI 18-40 and 40-60 keV or- Figure 3: Precessional 18-40/40-60 keV light 
bital light curves with hardness ratio (bottom plot). curves and hardness ratio. 

spectra are shown in Fig. |l| (right panel). It is seen that within errors they have an identical power- 
law shape with photon spectral index Fpi, ~ 3.8. Nevertheless, a tendency of the spectrum to get 
harder at ~ 0.25 (phase interval V in Fig. ^ middle panel) than at ^orh ~ 0-75 (phase interval I 
in Fig. H, middle panel) is clearly seen. Note also that at the middle eclipse (phase interval III in 
Fig. H, middle panel) the spectrum gets softer. This is exactly what is expected in the jet nutation 
picture: during the disk-jet nutation the angle between the line of sight and the jet axis changes 
by ~ 6 degrees, so the observer looks into the funnel less (at <porb = 0.5) or more (at <porb = 0.75) 
deeper, thus observing cooler or hotter parts of the jet base. 

3. Precessional variability. At a distance of 5.5 kpc to SS433, reliably derived from the kinematic 
properties of moving emission lines [||], the observed X-ray flux 18-60 keV corresponds to a max- 
imum uneclipsed hard X-ray luminosity of 3 x 10^^ erg/s. The precessional change of the X-ray 
flux is shown in Fig. ^. To plot this Figure, all available observations of SS433 by INTEGRAL 
with a total exposure time of about 8 Ms were used. This confirms the presence of a fairly broad 
region emitting in hard X-rays with size comparable to that of the accretion disk (r^ 10^^ cm), since 
(excluding orbital eclipses) the observed flux varies due to precession of the disk. The both preces- 
sional light curves show precessional variability and were used to constrain parameters of the hot 
corona ([^] and below). 

3. Results 

The INTEGRAL observations of SS433 provide three different light curves, which can be 
used to constrain the parameters of the system: 1) the orbital light curve, 2) the precessional light 
curve out of eclipses, and 3) precessional light curve in the middle of the eclipses (see Fig. ^, 
upper panels, and lower panels, blue and red crosses, respectively). We use a geometrical model 
described in detail in [^. Briefly, the model includes an optical star with mass My filling its Roche 
lobe, a compact star with mass Mx surrounded by an optically thick accretion disk with radius a^, 
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Figure 4: Joint analysis of orbital (upper plots) and precessional (lower plots) 18-60 keV light curves of 
SS433. On the bottom panels, precessional variability in the middle of eclipse (0.95 < (porb < 1-05) are 
shown by red crosses. 1). q — 0.1, 'short jet' corona (aj — 0.25, bj — 0.1, co — 40"). 2) q — 0.3, 'short jet' 
corona (aj = 0.35, bj = 0.13, « = 80"). 3) ^ = 0.3, 'short jet' corona {aj = 0.35, bj = 0.13, CO = 80"). 4) 
for the orbital hght curve for different mass ratios. 

and a hot corona which is modeled as a broad 'jet' parametrized by the part of ellipse with semi- 
axes aj and bj normalized to the binary orbital separation a. The elliptical corona is restricted by 
the cone with half-angle CO. The shape and amplitude of the precessional light curve constrain the 
height of the hot X-ray corona, while the orbital light curve restricts the accretion disk radius. At 
a given binary mass ratio, after finding the best parameters for the precessional variability, we can 
calculate the deviations of the model orbital light curve from the observed one. The results of the 
joint analysis of the orbital and precessional 18-60 keV light curves of SS433 are shown in Fig. ^ 
To avoid contamination from thermal X-ray emission from relativistic jets, we repeated the 
analysis using only hard X-ray 40-60 keV light curves shown in Fig. § and Fig. |3[ The result 
for different mass ratios q = M^/M^, is presented in Fig. It is seen that the hard X-ray orbital 
and precessional light curves can be simultaneously reproduced by the geometrical model for the 
binary mass ratio q > 0.3 (at smaller mass ratios a plateau corresponding to the total eclipse of the 
hot corona by the optical star appears in the orbital light curve, which is not observed), in agreement 
with the analysis of the broadband 18-60 keV orbital and precessional light curves shown in Fig. 

i 

4. Conclusions 

1) INTEGRAL observations of SS433 in hard X-rays 18-60 keV allowed us for the first time 
to make orbital-resolved spectroscopy of the X-ray eclipse in the precessional phase corresponding 
to the maximum opening angle of the disk. The hard X-ray continuum is fitted with a power-law 
with photon index F 3.8 which does not significantly change across the eclipse, suggesting the 
origin of this emission as being due to scattering in hot corona surrounding the funel around the 
jets in a supercritical accretion disk in SS433. 

2) For the first time, joint analysis of hard X-ray (40-60 keV) orbital and precessional light 
curves has been performed. This analysis independently confirms our previous result [||] that the 
low value of the mass ratio q = M^/My, in SS433 cannot reproduce the observed orbital and preces- 
sional light curves. With the existing estimates of the mass function of the compact star, the most 
likely value g ~ 0.3 points to the black hole nature of the compact star in SS433. 
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Figure 5: Joint analysis of orbital (upper plots) and precessional (lower plots) 40-60 keV light curves of 
SS433. Left panel: q = 0.1, 'long jet' corona (aj — 0.25, bj — 0.55, CO = 80"); only orbital light curve can 
be reproduced. Middle panel: ^ = 0.1, 'short jet' corona (aj = 0.25, bj — 0.1, co = 80"). Both orbital and 
precessional light curvs can be fitted, but the total eclipse (plateau at zero flux) appears (unobserved). Right 
panel: q = 3, 'short jet' corona (aj = 0.35, bj = 0.13, CO = 80"); both orbital and precessional light curves 
are well reproduced. 

3) The shape of the hard X-ray orbital light curve 40-60 keV demonstrates two humps at 
around orbital phases 0.25 and 0.75, likely due to the nutation effects in SS433. 

New INTEGRAL observations of SS433 at different precessional phases will be used to further 
constrain physical parameters of this unique galactic microquasar. 
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